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We investigate entanglement of spin pairs in alternating open spin chains in the equilibrium state 
in an external magnetic field. We calculate the reduced density matrix of spin pairs and estimate 
the concurrence with Wootter's criteria. The obtained results demonstrate the dependence of the 
entanglement on the temperature, chain's length, the positions of the spins, and the ratio of the 
spin-spin interaction constants. 
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I. INTRODUCTION 

A study of entangled states is an important direction 
in quantum information theory. A quantitative investi- 
gation of the degree of entanglement often poses consid- 
erable difficulties and is possible only in relatively simple 
cases. The problem becomes much simpler, if it is known 
how to diagonalize the Hamiltonian [1-3]. In particular, 
many interesting results were obtained for spin rings and 
infinite spin chains [4-6]. On the other hand, possible 
experimental realizations of qubits [7] involve finite open 
chains which were also used for the transfer of quantum 
states from one end of the chain to the other [8, 9]. Im- 
portant results were obtained by Wang [10] who investi- 
gated open boundary effects on the ground-state entan- 
glement in a homogeneous Heisenberg model. However, 
the models [4-6, 10] do not allow one to solve the qubit 
addressing problem which can be attacked in systems 
with several different Larmor frequencies. Such Larmor 
frequency variations naturally emerge in inhomogeneous 
systems with different constants of the spin-spin interac- 
tion (SSI) for different pairs of neighboring spins. 

In this paper, we investigate entanglement in open al- 
ternating spin chains. Open alternating spin chains with 
XY-Hamiltonians can be used as a first step for the so- 
lution of the qubit addressing problem. Recently, the 
XY-Hamiltonian of an open finite alternating spin chain 
has been diagonalized [9, 11] on the basis of a generaliza- 
tion of the classical methods [1] for homogeneous chains. 
It was shown that the transfer of quantum states from 
one end to the other is possible with 100% fidelity for 
longer alternating chains than in the case of homogeneous 
chains [9]. 

The concurrence [12] is used as a quantitative measure 
of entanglement. A possibility of the exact diagonaliza- 
tion of the Hamiltonian [9] makes it possible to construct 
the reduced density matrix for any spin pair in a chain 
described by the equilibrium density matrix. Then entan- 
glement can be investigated with Wootter's criteria [12]. 
As a result, one can study the influence of the temper- 
ature, the chain length, the distance from the ends, and 
the ratio of the SSI constants on the entanglement of any 
spin pair in the chain. 

The present paper is organized as follows. In Section II 



the model Hamiltonian of the open alternating spin chain 
is described and the eigenvalues of the Hamiltonian and 
its eigenvectors are given. In Section III we obtain the re- 
duced density matrix for all pairs of the neighboring spins 
of the chain. We study the dependence of the spin pair 
entanglement on the temperature, chain's length, spin 
positions and ratio of the SSI constants in Section IV. 



II. THE MODEL 

We consider an open alternating spin-1/2 chain with 
the XY-Hamiltonian in a strong external magnetic field. 
The Hamiltonian of the system can be written as 



N 

H = ^ u n I nz 

N-l 

^ Dri,n+1 (-^n.x-^n+l.x ~t~ -^n,y-^n+l,y)- (1) 



Here, I na are the spin-1/2 matrices (a — x,y,z), and 
N is the number of the spins. The alternating Larmor 
frequencies are determined by the equation 



u>\, n is odd 
u>2, n is even. 



(2) 



The alternating coupling constants are 



D 



n,n+l 



Di, 
D 2l 



n is odd 



n is even. 



(3) 



We assume that the number of the spins, N, is odd. The 
Hamiltonian of Eq. (1) can be diagonalized [9, 11] on the 
basis of a generalization of the classical methods [1] for 
homogeneous chains. 

The Jordan- Wigner transformation [1, 2] that maps 



2 



spins onto free spinless fermions can be written as 

= In, x - U n ,y = (-2)"- 1 ([] I, A Cn, (4) 



\l=l 
(n-\ 



In,+ = In, X + U n ,y = (-2)"" 1 ]J hz 4, (5) 



In,z — c n c n 1/2 



(6) 



where and c„ are creation and annihilation operators. 
As a result, the Hamiltonian in the matrix notation has 
the following form 



1 1 N 

H=-c\D + 2Q)c--J2LO n 



(7) 



n=l 



Here we denote the row vector — (c\, . . . , cjy), the 
column vector c = (ci, . . . , cjv) (the superscript t repre- 
sents the transposition) and D + 2f2 is a three-diagonal 
matrix. 

The diagonalization of the matrix D + 2f2 is performed 
by the unitary transformation 

D + 2il = UA U\ A = diag{Ai,...,Ajv}. (8) 

The new fermion operators 7^ and 7^, which are intro- 
duced by the relationships 



N 



N 



(9) 



fc=i 



k=l 



transform the Hamiltonian of Eq. (7) into the Hamilto- 



nian 



1 N 1 N 



(10) 



fc=l 



n=l 



This is a free fermion Hamiltonian with the energy lev- 
els 1/2 A/j. Finally the eigenvalues X k and eigenvectors 
\uk) = (uik,U2k, ■ ■ ■ , UNk) 1 of the Hamiltonian of Eq. (10) 
can be written as [11] 



'w! + uj 2 + ^/(ujx -uj 2 ) 2 + D\A kl 



Afc 



where 



1,2, 



jV-1 

2 



2wi, fc = ^ 



(11) 



jV+3 jV+5 



,N 



A k = l + 25cos(^^j+5 2 , 5 = D 2 /D 1 . 



(12) 



For all the indices fc = 1, . . . , N except the index k = 



(N + l)/2 the eigenvector \u k ) has the elements 



u 3,k = < 



4, Jl 

^ A fc -2 Wl 



5 sin ( ZL 7 ^p^ 



j = l,3,5,...,N 
A fe sin(^), 
U = 2,4,...,JV-1 

with the normalization coefficient 
2|A fc -2wi| 1 



(13) 



VW+T V(Afc-2 Wl ) 2 + r>?A fc " 

The elements of the eigenvector |tt(jv+i)/2) rea d 



(14) 



%,(JV+l)/2 

with the normalization coefficient 



i?M)(^)/ 2 , i = 1,3, 5,..., N 
0, .7 = 2,4,. ...JV-l 



5 



1 



5N+1 _ 1 



1/2 



(15) 



(16) 



For D\ = D'2 = D and cji — u>2 — loq we obtain a 
homogeneous chain and the Hamiltonian of Eq. (10) can 
be written as [3] 



N 1 
H = H X kl\lk - ^Nujq, 



(17) 



fc=i 



with Afe = Dcos + u>o and eigenvector \u k ) has the 
elements 



Uj.k 



N+l 



1/2 



sin 



irkj 
N+l 



.7 = 1,2,..., N. 



(18) 

For homogeneous chains these expressions are valid both 
for odd and even N. 



III. THE REDUCED DENSITY MATRIX 

In this section we consider a one-dimensional many- 
spin system in a thermodynamic equilibrium state. We 
describe an algorithm for obtaining the reduced density 
matrix for any spin pair. This algorithm will be ap- 
plied in order to obtain the reduced density matrix of 
the nearest-neighbor spins i and i + 1 for the open alter- 
nating chain with an odd number of the spins. 

The density matrix, p, in the thermodynamic equilib- 
rium system is 



-(3H 



P = 



(19) 



3 



where (3 = h/kT, T is the temperature and Z = 
Tr{e~@ H } is the partition function. The density matrix, 
p, can be written as 



p= E 

?l,«2,...,«f 



J N ' 



(20) 



where N is a number of spins in the system, (k = 
1, 2, . . . , N) is one of the values {0, 1, 2, 3}, x° k = I k is the 



unit matrix of the dimension 2 x 2, x^ = Ifc x , x 2 , = 7, 



ky ■ 



x k = Ikz, and ck*2 jv * s a numerical coefficient. 

In order to obtain the reduced density matrix for spins 
we consider the system of all other spins as the envi- 
ronment. Averaging the density matrix, p, over the envi- 
ronment and taking into account that Tr{a;|' c } = (k = 
1, 2, . . . , N; £fc = 1, 2, 3) we find for the reduced density 
matrix, p, of spins i and j the following expression 



Pij 



3 

E 



.7:7 



(21) 



where the coefficient af/ 3 is defined as 



o£ 4j = 



2 Ar - 2 Tr{/ M ,f a A} 
Tr{(^)2(4) 2 } ' 



(22) 



The calculations of the traces in Eq. (22) are performed 
in an 2 Ar -dimension space. Here, the density matrix, p, 
and the density matrix, p^ , are normalized to unity. 
Since [77, I z ] = one can see that 



a 01 = a 10 



a 



02 



i 20 -- 

^13 



0. 



,23 



,,31 



= <4f=0. (23) 



.12 - 



= 4Tr {pl ix l jy } 



0, af = 



Furthermore, 
4Tr {pl iy l jx } 
4Tr {phyljx} = a- 

As a result, the structure of the reduced density matrix 
of two spins is given by 



= 0, and a# = 4Tr{p7 ix 7 jx } = 
J2 from the symmetry of the problem. 



Pij 



fa 0\ 
6x0 
x c 

\0 



(24) 



wherea= 1/4 + ^ + ^ + ^,6= 1/4-^ + ^-^, 

„03 30 33 03 30 33 

c= 1/4 + ^-^-^=1/4-^-^ + -^ and 

33 and 



x = 



:so 



— -7?-. We calculate the coefficients ex 03 , a;-.- 
a}} below. 



a) 



An open alternating chain with an odd number 
of spins 



In order to calculate the coefficients a 



03 



v30 



„33 



and 



alj we use the diagonal form of the fermion Hamiltonian 



of Eq. (10) of the system. One finds that 
< = 1/4, 



a 3 ° = Tr{p/ i2 } = ^ U 2 p - 



a% = Tr {pl jx } = 5>? p - 



e 2 



+ e 2 a p 



1 

2' 



e 2 



+ e 2 a p 



A P 2' 



(25) 
(26) 

(27) 



a 33 =4Tr{p7^7 jz } 

= 4Tr{p (r.r.r r, - icjc; - ^cjcj + |) } 




e 2 - 



1 + e" 




- 2 £«; 



1 + e -2 A *,- 



e 2^p 



+ e 2 a p 



(28) 



Our numerical simulations for the spin chains with 9 spins 
demonstrate that entanglement appears only between the 
nearest-neighbors even if SSI of all spins are taken into 
account. The following coefficients are sufficient for our 
investigations: 



n 

*»,*+! 



22 



+1 = 4Tr {pli x li+i, x } 



Tr {p( c 14+i + 4^+i 



+i)} 



— 2 ^ ' Ujn u i+l,n 
n 1 

Introducing the notations 



e 2- 



+ e 2- 



(29) 



Ci 



w 2 — Wi 
_ Di - : 



C 2 



W2 + Wl 

_ Di - ' 



Lk 



1 7. - 1 JV=i 



I u _ N+3 

' — 2 ' ' ' ' ' 



C-L-y/Cf+Ak 



N 



(30) 



fk = A\ 



N+l 



(C 1 + y/cf+A^) 2 +A h 



1,... 



1 



JV-1 



W+:L I " (Ci-VC 1 2 +A fc )2 + A fc 
3 

2 



(31) 



4 



R k = Ai 



2 / D 1 



fc — -L, . . . , 2 



L _ W+3 at 
— 2 , . . . , JV 



(32) 



C 2 + y /Cl+A k , k = 1, . . . , 

C 2 - v / Ci 2 + A fc J fc = ^,...,iV 



and 



5(efe) = 



1 + e -1 "^ 



(34) 



we obtain explicit expressions for the coefficient oif\ +1 . 
Since the alternating chain is non-symmetric we obtain 
the different expressions for the coefficients at even and 
(33) odd i. At even i the coefficient is 



sin 



■Kvni 

FT 



N + 1 J \ Di 



+ [ J2f n L n Ss(n) sin (^q^J J XI f m L m S s (m) sin ^-^-Jg(e m ) 



r. 



7TTTL1 \ . ( TTUl , 

- I sin I - . : |ff(em) 3(e«) 



- X X fmfnL n L m Ss(m}S s (n) sin . ^ 1 , _ i jy , x 

/ N + l \ / .. / JV + 1 \ 1 



r 



+ 1, (35) 



where S s (m) = [,5 sin + sin (^±21) 



At odd i the coefficient af^ +1 is 



<*M+i = 4 



^ ^ f n R m sin 2 (^)«Wff(gj(e,)+ E /-^("f^W f^ 11 )* (^) 



^/nin<55(n)sin f-^-pY^j ] ^2fmL m Qs(m)sm y ^+i jdi^ 



E E fmfnLnL m Qs(m)Q s (n) sin ( - 

l N + l n-Am \ 



7rm(i + l)\ . ( im(i + 1) 



7V+ 1 



2<j ^ M^)+B 2 Hr'j(y)+ E / P sin 2 (^^) 5(ep )Ul, (36) 



-hi | \g(e m ) g{e n ) 

. 2 / + 1) 



where Q s {m) = [<5sin ( 7r ^ 1) ) + sin . Ex- can be obtained in a similar way. 

plicit expressions for the coefficients a° 3 i+1 , af° +1 and 



■5 



B. The open homogeneous chain 

The expressions of the previous Section simplify for 
homogeneous chains and can be written as follows 



„03 



a 30 



N + 



sin 



N + 



sin 



ink \ . . 1 ,„„. 
T JS-(efc) - ^ ( 38 ) 



7V + 



cry = 27 y . The concurrence of the two-spin system with 
the density matrix pab is [12] 

Cab = max{0, 2A — Ai — A2 — A3 — A4} 

A = max{Xi, A 2 , A3, A4} (42) 

where Ai, A2, A3, and A4 are the square roots of the 
eigenvalues of the product pabPab- Since both pab and 
Pab are positive operators, it follows that their product, 
though non-Hermitian, has only real and non-negative 
eigenvalues [13]. 

Using Eq. (24) one can find the expressions of Ai, A2, 
A3, A4 as follows 



a 33 = 



16 



(N + l) 



sEE shr 



k p^k 

■ 2/' j^P 



sin 



N+l 



ink 
N+l 



Xi = A 4 = Vad, 



A 2 , 3 = \x±Vbc\. (43) 



EE 

k p^k 



ink . inp 
sin sin ■ 



7V + 1 N + l 



. jnk . jnp \ 
sm __ sm __ 5 ( efc ) 5 ( ep) | 




sin 



ink 
N+l 



g(ek) 



+ E- 2 (] v tt).9M) + ^ (39) 



4 • / ink \ ■ f(i + l)Trk\ 

- sin ! — — - i Mi) 



k 



(40) 



Expressions for Eq. (43) yield all information necessary 
for the analysis of the pairwise entanglement for such 
chains. 



V. NUMERICAL ANALYSIS OF THE 
CONCURRENCE IN SPIN PAIRS 

As mentioned above, our numerical calculations in the 
nine-spin chain with the Hamiltonian of Eq. (1) lead to 
the conclusion that the concurrence is non-zero for near- 
est neighbors only We focus below on entanglement of 
the nearest neighbors in the open alternating chain with 
the XY-Hamiltonian in the case of zero Larmor frequen- 
cies. 



where D\ = D 2 = D, and the one-fermion spectrum 
is defined as 

e, = 2cos^ T + ^, (fc = l,2,...,7V). 

We emphasize again that Eqs. (37)-(40) are valid for 
all N. 



IV. MEASURE OF ENTANGLEMENT 

In this paper we restrict our attention to the analysis 
of the entanglement between two arbitrary spins in the 
chain. Here we consider the concurrence as a measure of 
entanglement [12]. 

Let A and B be a pair of qubits, and let the density 
matrix of the pair be pab which may be pure or mixed. 
Then the "spin-flipped" density matrix is determined as 




PD/2 



PAB = K O O y )p*A B (O y ® (Ty) 



(41) 



where the asterisk denotes complex conjugation in the 
standard basis {|00), |01), |10), |ll)}and the Pauli matrix 



FIG. 1: The concurrence of spins 2 and 3 versus the temper- 
ature for the number of spins N — 101 and D2/D1 — 1.5. 
The inset displays the concurrence versus the temperature on 
a greater scale. 

Using Eqs. (25)-(40), (43) we investigate numerically 
the dependence of the concurrence, C12, of the first and 
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FIG. 2: The nearest-neighbor concurrence versus the site 
number in the homogeneous chain for N = 100. The concur- 
rence oscillates with a two-site period. The oscillations decay 
when the spin pair is far from the ends of the homogeneous 
chain. 



second spins on the temperature. Figure 1 shows that 
entanglement appears at (3D\ w 1. This corresponds 
to the temperature T rts 0.5 fiK at D\ w 27r • 10 4 c _1 . 
The temperature at which entanglement emerges in the 
pairs of neighboring spins depends on the ratio of the SSI 
constants, chain's length and the distance from ends of 
the chain. It is interesting to note that ordered states 
of nuclear spins were observed at microkelvin tempera- 
tures [14]. 

Numerical results show (Fig. 2) that entanglement os- 
cillates with the two-site period. The qualitative expla- 
nation of the numerical results is the following. Qubits 1 
and N are situated at the chain ends. Since entanglement 
is non-zero only between the nearest-neighbors, the pair 
of qubits 1 and 2 and the pair of qubits N-l and N have 
the maximal pairwise entanglement for the homogeneous 
chain. Spin 2 can be entangled both with spin 1 and with 
spin 3. Since spin 2 is strongly entangled with spin 1, the 
entanglement of spins 2 and 3 is weaker. As a result, spin 
3 is strongly entangled with spin 4, etc. This explains the 
oscillator behavior of the concurrence displayed in Fig. 2. 
The oscillations decay when the spin pair is far from the 
ends of the homogeneous chain. These oscillations do 
not decay for alternating open chains as it is displayed 
in Fig. 3 at different ratios of the coupling constants, 
DijD\. The oscillation of the entanglement from zero 
to values close to one is due to different coupling con- 
stants between spins in the alternating chain. In fact we 
have here the dimerised spin chain which can be consid- 
ered qualitatively as a set of non-interacting spin pairs at 
D2/D1 > 2. The dependence of the concurrence on the 
ratio of the SSI is shown in Fig. 4. It is worth to notice 
that it is impossible to observe the influence of ends of 
the chain on the entanglement already at D2/D1 = 1.5. 




FIG. 3: The nearest-neighbor concurrence versus the site 
number for N — 17 and f3D\/2 — 30 at the difference ratios of 
the SSI coupling constants; & :D2/D\ = 1; h-.DijDi = 1.17; 
c:D 2 /Di = 3. 
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In Figs. 5-6 the dependence of the concurrences C12 
and C23 on the length of the homogeneous chain consist- 
ing of an odd (even) number of spins is displayed. The 
concurrence C\i decreases at even N and increases at odd 
N when N increases. For the concurrence C23 the situ- 
ation is opposite. In order to explain these results it is 
necessary to take into account the influence of the ends of 
the chain on the entangled states. As a result, we obtain 




FIG. 4: The nearest-neighbor concurrence of the pair of spins 
1 and 2 and of the pair of spins 2 and 3 versus the ratios of 
the SSI for N = 55 and for 0Di/2 = 30. 



that the ends of the chain yield opposite contributions 
to the entanglement for the chain of an odd number of 
spins. For example, one can easily find in a chain con- 
sisting of 5 spins that spin 5 leads to a decrease of the 
entanglement of spins 1 and 2. This effect is diminished 
when the odd number of spins increases and then the en- 
tanglement of spins 1 and 2 increases. On the contrary, 
the second end of the chain increases the entanglement 
of spins 1 and 2 at an even number of spins. This leads 
to the calculated decrease of the concurrence (see Fig. 5) 
when the even number of spins increases. Analogously, it 
is possible to explain the growth (fall) of the concurrence 
of the entanglement of spins 2 and 3 in the dependence 
on the length of the chain consisting of an even (odd) 
number of spins (see Fig. 6). The maximums of the con- 
currence in Fig. 6 have a simple qualitative explanation. 
The growth of the concurrence up to its maximal value 
is determined by the decrease of the partial suppression 
of the SSI due to the flip-flop transitions at low tem- 
peratures. The state of spins 2 and 3 approaches to a 
separable one at low temperatures when the spins are 
trying to orient parallel to the local fields. Entanglement 
of such a state decreases when the temperature decreases. 





D 2 /D=1 



10 

(3D/2 



3 spins 

7 spins 

4 spins 

6 spins 

o 100 spins 

+ 105 spins 

20 25 



FIG. 5: The concurrence of spins 1 and 2 versus the tem- 
perature for the number of spins N — 3,4, 6, 7, 100, 105. The 
behavior of the concurrence C12 is different for chains with 
even and odd N. 



3 spins 

7 spins 

4 spins 

6 spins 
o 100 spins 
+ 105 spins [ 




15 

PD./2 



FIG. 6: The concurrence of spins 2 and 3 versus the temper- 
ature for the number of spins N = 3, 4, 6, 7, 100, 105. 



VI. CONCLUSION 

In the present work we investigate entanglement in al- 
ternating open chains. The developed methods of the di- 
agonalization of the XY-Hamiltonian of open alternating 
spin chains allow us to study the pairwise entanglement 
for different parameters of the chain and its temperature. 

Similar methods can be applied for an investigation of 
the pairwise entanglement in open chains with periodic 
coupling constants [15]. Methods of exact diagonaliza- 
tion of the XY-Hamiltonian in such chains [15] allow a 
solution of different problems of quantum information 
theory for models of quantum registers taking into ac- 
count the qubit addressing. However, the qubit address- 
ing requires a large difference of the Larmor frequencies of 
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different spins in comparison with their SSI coupling con- 
stants. Our preliminary calculations demonstrate that 
the concurrence is close to zero at such conditions and 
the entangled states do not appear. The question is a 
subject of our further investigations. 
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